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ABSTRACT 

An understanding of the behavior of vegetation canopy reflectance as a function 
of solar zenith angle is important to several remote sensing applications. Spectral 
hemispherical-conical reflectances of a nadir looking sensor were taken throughout 
the day of a lodgepole pine and two grass canopies. Mathematical simulations of 
both spectral hemispherical-conical and bi-hemispherical reflectances were per- 
formed for two theoretical canopies of contrasting geometric structure. These re- 
sults and comparisons with literature studies showed a great amount of variability 
of vegetation canopy reflectances as a function of solar zenith angle. Explanations 
for this variability are discussed and recommendations for further measurements 
are proposed. 


*This is a preprint of an article submitted to Photogrammetric Engineering and Remote Sensing. 
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INTERPRETING VEGETATION REFLECTANCE MEASUREMENTS 
AS A FUNCTION OF SOLAR ZENITH ANGLE* 


INTRODUCTION 

An understanding of solar radiation interaction with vegetation canopies is necessary for ac- 
curately interpretating remotely sensed data. Broad and narrow-band spectral reflectance measure- 
ments of vegetation canopies are often used to characterize this solar interaction. Measured re- 
flectance however, is a complex function of canopy constituent optical properties (Gates, 1970; 
and Knipling, 1970), canopy geometry (Kimes, ct al., 1979b; and Ross, 1976), optical properties 
of the ground, atmospheric conditions (Kriebel, 1976; and Ross, 1976), solar zenith angle (Smith, 
et al., 1974; Kriebel, 1975; and Jarvis, et al., 1976), and sensor inclination and azimuthal view 
angles (Smith and Oliver, 1974; Kriebel, 1978; and Smith, et al., 1979). 

Understanding canopy reflectance as a function of solar zenith angle is important for several 
remote sensing applications. For example, such knowledge can improve multitemporaj vegetation 
classification by using sun-angle signature extension techniques (Smith, et al., 1975). At higher 
latitudes low sun-angles predominate and an understanding of the reflectance changes at low sun- 
angles would be beneficial. Diurnal reflectance trends are also important in photosynthetic and 
productivity studies. 

To better understand these relationships, spectral reflectance measurements were obtained for 
several solar zenith angles for a lodgepolc pine and two grass canopies. Mathematical simulations 
of the diurnal reflectance from theoretical vegetation canopies were performed. The instruments 
and methods used to obtain these data are described and the resulting trends presented and dis- 
cussed. The results were also compared to other field measurements of vegetation canopies, and 
sources of variation are discussed. 


*This is a preprint of an article submitted to Photogrammetric Engineering and Remote Sensing. 
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REFLECTANCE NOMENCLATURE 

The two reflectance measurements which arc most commonly reported in the literature for 
natural vegetation canopies are bi-hemispherical and hemispherical-conical reflectance. In this 
study a nadir looking sensor was used to measure the hemispherical-conical reflectance of vege- 
tation canopies. The definitions of these reflectances are presented as follows. For clarity the 
standard nomenclature and symbolism for the basic radiometric quantities (c.g., radiance, irradiance, 
and exitance) as presented by Suits (1975) were used exclusively. 


Bi-hemispherical reflectance (p 11 ) is defined as the ratio of the reflected exitance to the ir- 
radiance at the target surface. The hemispherical-conical reflectance (p c ) for a nadir looking 
sensor having a field of view of less than 27r steradians is measured as the ratio of the reflected exi- 
tance of a surface in the direction of the sensor’s field of view (FOV) to the reflected exitance of a 
perfect .reflecting horizontal lambertian surface in the direction of the sensor’s FOV. The mathe- 
matical representation of the above measurements is given by: 
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where 0 r , 0 r = The zenith and azimuth angles of the reflected sources respectively 
0 j, <j>\ = The zenith and azimuth angles of the incident sources, respectively 

Lp Lj ® The radiance values of the reflected and incident sources respectively as a 
function of 0 and <j> 
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where A = The denominator of Equation (1). 

0, ,02,01 ,#2 “ The azimuth and zenith angle limits of the sensor’s steradian view. 

Note that the denominator and numerator of Equations (1) and (2) represent what tiie sensor 
would measure from a reference panel (e.g., barium sulfate) and vegetation canopy respectively. 


The bi-hemisphcrical (p 11 ) and hemispherical-conical (p c ) measurements are related to the 
bi-directional reflectance function (p B ) as described by Nicodemus (1970) and Kriebel (1976). The 
p B function is 


P u (0p ^r) 


d L r (0 r , 0j ■) 
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(3) 


As stated by Kriebel (1976) the function is defined as the relation of that part of the total spectral 
radiance dL r (0 r , 0 r ) reflected into the direction 0 r , 0 r which originates from the direction of in- 
cidence 0i,0i, to the total spectral irradiance LjfOi^) cosOidWj impinging on a surface from the 
direction 0;,0 5 . This particular bi-directional reflectance function (p B ) is a unique characterization 
of a surface and is not dependent on the irradiance distribution as are a number of other bi- 
directional functions presented in the literature. 


From Equation (3) it follows that 
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As a consequence, p c and p 11 are dependent on both the bi-directional reflectance function p B 
and the solar irradiance distribution (Equations (1) and (2)). From these equations it is important 
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to note that hemispherical-conical reflectance (p c ) which is most commonly utilized in remote 
sensing research is not a physical parameter uniquely determined by a particular vegetation canopy 
but is dependent on the anistropic distribution of sky irradiancc. it is the bi-directional reflec- 
tance function together with information on the anistropic distribution of sky irradiancc which 
ultimately determines how vegetation canopy p c and p 1! reflectances will behave under various 
solar irradiancc conditions. Bi-directional reflectance is a function of a number of geometric and 
optical characteristics of the canopy components as discussed by Oliver and Smith (1974) and 
Ross (1976), Kriebcl (1974, 1978) presents spectral bi-directional reflectance values for four 
vegetation canopies. 

A description of the above and other types of reflectance measurements are presented by 
Judd (1967) and Nicodcmus (1970). The reflectance measurements can be distinquished as broad 
band (p H ,p c ,p li ) or spectral (p” , p^.pjj) reflectance. The variability within and between these 
reflectance measurements are explored below. In addition, throughout this paper p£ will refer to 
a nadir sensor angle. 

INSTRUMENTATION AND METHODS 

Both field and simulated data were obtained. In all field experiments a scene recording 
radiometer (SRR) as described by Berry, et at., (1978) was used to obtain spectral hemispherical- 
conical measurements (p^) of vegetation canopies. The SRR was suspended on support cables 
attached to two !5in towers which allowed nadir looking measurements with a field of view of 
22.5° from above the canopy to be obtained. The optics consisted of a six narrow band inter- 
ference filter wheel Interfaced to a Hasselblad EL500 camera which provided a photographic 

record of the scene. All filtered spectra! data were referenced to a horizontal barium sulfate panel 

\ 

to provide reflectance values. Filters used were centered at 4800, 5500, 6750, 7300, 8000, and 
9600 A and had a half width bandpass of 100 A. 
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Four experiments were performed at two field sites to evaluate the spectral hemispherical- 
conical reflectance (p£) of various targets with changing solar zenith angle. The targets and date 
of measurements at site I were: lodgepolc pine with a grass understory (August 4, 1976) and open 
meadow (August 6, 1976). At site 2 a grass community was measured (April 20, 1978 and May 
18, 1978). In all experiments the sensor was nadir looking and a constant sensor position on the 
tower system was maintained. Target radiance and barium sulfate panel radiance data were 
taken frequently throughout the day and the corresponding measurements were ratioed to obtain 
the values. 

A description of the two study sites was as follows: Site 1 was located southeast of Lcad- 
vilie, Colorado, in the northeastern section of the Iron Hill area at an elevation of approximately 
3,322m. The 45.7m long transect defined by the tower system was oriented at 60° and main- 
tained a constant slope and aspect of 5% at 045°. The vegetation gradated from a relatively dense 
lodgepolc pine stand to an open, grass-covered clearing. Average height of the iodgepole pine 
stand was approximately 6.1 m. Canopy density was variable, ranging from 80% crown closure to 
the open meadow previously mentioned. The meadow was populated with rush ( J uncus sp. ) and 
sedge ( Carcx sp. ). A more detailed description of the site was presented by Heimes and Smith 
(1977). A description of the geometric structure of the lodgepolc pine canopy was presented by 
kimes, et al., (1979a), 

Site 2 was located west of Fort Collins, Colorado, at an elevation of 1,570m at the Colorado 
State Forest Service Nursery. The towers at each end of the transect were 9m high. Vegetation 
along the transect consisted of grasses which included fescue ( Festuca ), bluegrass (Poa), sedge 
(Carex), wheatgrass ( Agropyron ) and brome (Bromus ). A more detailed description of the site 
was presented by Ranson, et al., (1978). 

A solar radiation vegetation canopy (SRVC) model was utilized to simulate the p£ and p” 
reflectances for various theoretical vegetation canopies. The SRVC model is a Monte Carlo model 
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which physically accounts for variations in direct/diffuse solar irradiance ratios, r.oiar zenith angle, 
leaf angle geometry, leaf area indices (LAI), leaf spatial dispersion, and leaf and soil optical proper- 
ties (Oliver and Smith, 1974). Several vegetation canopies were simulated with various geometric 
structures and LAI for a photosynthetically active spectral wavelength (0.68pm). The entire pro- 
cedure and canopy absorption results are presented by Kimcs, ct ai., (1979c). Two canopy geom- 
etries were simulated: an erectophilc (mostly erect leaves) and planophilc (mostly prostrate leaves) 
as described by DcWit (1965). These canopies were simulated for LAI of 1,0, 4.0, and 7.0. The 
simulated results were compared with field measurements and arc presented in this study. 

Finally, literature studies on vegetation reflectance trends as a function of solar zenith angles 
were compared with the results of the above experiments. 

RESULTS AND DISCUSSION 

Figures la and b present the 0.68 and 0.80pm band hemispherical-conical reflectances (p£) 
as a function of solar zenith angle for the lodgepolc pine and meadow at site 1. These data were 
acquired on August 4, 1976, from 0657 to 1400 hours mountain daylight time (MDT) and August 
6, 1976, from 0815 to 1558 hours MDT, respectively. The arrows denote the sequence of data 
points from morning to afternoon. For both canopies p^ increased with decreasing solar zenith 
angle except for the 0.80pm band of the meadow in the morning. The meadow values were 
lower in the afternoon with respect to tire morning and the opposite trend was apparent for lodge- 
pole pine. AI! of the above p^ measurements at site i were taken when the direct solar path was 
free of clouds. Consistently, there was a build up of cumulus clouds around noon which continued 
into the afternoon. The experiments Were terminated when cloud cover did not permit a direct 
solar path free of clouds. 

Figure 2a and b presents the measured p£ results as a function of solar zenith angle for the 
grass community at site 2 on April 20, 1978, and May 18, 1978. The April and May data were 
collected from 0856-1715 MDT and 1015-1947 MDT, respectively. Data for both the 0.68pm 
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und 0.80pm bands arc presented. Both of these experiments were executed under rapidly changing 
sky conditions. Cumulus cloth)? partially or totally obscured the sun at various times throughout 
the measurement periods. Accordingly, there was a relatively large amount of variability in the 
data. No dominating trends occurred as a function of solar zenith angle with exception of the 
0.68pm band on May 18, which showed a slight increase in reflectance with decreasing solar zenith 
angle. 

The SRVC simulated p£ reflectance for a nadir looking sensor and the p” reflectance for an 
crcetophilc and planophile canopy of three LAI values are shown in Figure 3a, b, c, and d. 

These results showed that for LAI of 1.0 and 4.0 of the ercctophilc canopy the p^ decreased with 
increasing solar zenith angle (Figure 3a). However, the p” for all three LAI values tended to in- 
crease with increasing solar zenith angle (Figure 3b). These differences arose from anistropic scatter- 
ing by the vegetation canopy into off-nadir view angles. The effect can be correctly quantified by 
the bi-directional reflectance function as presented by Kriebcl (1976, 1978) and discussed below. 

In contrast, the p^ for the ercctophilc canopy of an 7.0 LAI increased with increasing solar zenith 
angle (Figure 3a). Subtle geometric and optical effects of the soil and canopy accounted for this 
discrepancy as will be discussed later. 

Grass canopies such as measured in this study tend to assume crectophiic geometries (Oliver 
and Smith, 1974). The measured results of the meadow (Figure lb) showed p^ both increasing and 
decreasing with increasing solar zenith angle. 

The geometry of todgepole pine was intermediate between a planophile and ercctophile canopy 
as measured by Kirnes, et a!., (1979a). Kimcs, et al., (1979b) have used the SRVC model to simu- 
late p£ as a function of solar zenith angle for the specific geometric and optical properties of the 
lodgepole pine canopy at site 1. Simulated p£ decreased slightly as , v.olar zenith angle increased 
as was shown by the measured data (Figure la). 
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Tlie simulated results of the planophiie canopies were much less variable than the erectophilc 
canopy (Figure 3c, d). 

The diurnal reflectance results from this study were compared to previous studies. These 
comparisons showed a great deal of variability in diurnal reflectance trends. This variability and 
explanations for it will be presented under the headings of p c , p" , and p 11 reflectance. 

(f Reflectance 

In general the measured and simulated p£ results for the vegetation canopies supported a 
trend of decreasing p£ with increasing solar zenith angle. There were exceptions, however, as noted 
above. 

Smith, et ul., (1975) have measured p£ of wheat "Using an Exotech ERTS Radiometer. The 
spectral rcfler .nyt A three plots were sampled using the four Landsat MSS bands. Mean p£ 
values for various solar zenith angles were recorded for three different wheat stages: tiller- 
ing (March 20), jointing (April 23), and heading (May 20). Figure 4a, b, and c showed the mean 
p£ values of each wheat stage for bands 5 and 7 (0.6-0.7pm and 0.8-1. 1pm, respectively) as a 
function of solar zenith angle. In addition, the mean LAI of each stage was presented. Although 
the data in Figure 4a, b, and c were collected for only a few solar zenith angles, the results showed 
temporal variation of the diurnal trajectories of wheat. The measured geometric structure of the 
three stages of wheat approximated erectophilc geometries (Smith, ct al., 1975). The diurnal tra- 
jectories both increase and decrease with increasing solar zenith angle. Explanations for this vari- 
ability are as follows. 

Kriebel (1978) has shown that the bi-directional reflectance functions for vegetation canopies 

were generally non-isotropic in nature, As a consequence, p^ variations arise from a changing 

anistropic irradiance field since the observed p£ is an integration of the bi-directional reflectance 
values (p^) for the near nadir reflectance angles and all incident radiation angles (Equations 2 

and 4). Each of these p)J values must be weighted by the distribution of the incoming solar 
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irradiance field. This weighting phenomena may account for the wide variability in measured 
seen on partly cloudy days (Fig. 2a, b). It is precisely this effect that is responsible for the chang- 
ing as a function of solar zenith angle on clear days. Specifically, the data from Kricbcl (1978) 
indicate non-lambertian behavior for the four vegetated surfaces: savanah, bog, pasture land, and 
coniferous forest. Assuming that the sun is the primary radiant source on a clear day, the p£ 
can be approximated by pj[(0 s ,0 s ;0 o ,0 o ) where 0 s> <j> s denote the solar zenith and azimuth angles 
respectively and 0 O and <p 0 denote the viewing angle of a nadir sensor (Equation 3). The 
pj[(0 s ,0 s ;0 o ,0 o ) values for the four vegetation types did not follow any consistent trend with 
increasing solar zenith angle. Thus from the results of this study and the cited literature one 
may expect p£ trajectories of vegetation as a function of solar zenith angle to be relatively vari- 
able due to the different bi-directional functions of plant canopies and variation in atmospheric 
conditions throughout the day. 

There was a significant deviation between morning and afternoon reflectances. In some in- 
stances the morning were higher than those in the afternoon (meadow, Figure lb) and in other 
instances the opposite is true (lodgepote, Figure la). Such an effect was not explored using the 
SVRC model because azimuthal symmetry is mathematically assumed, Ripley and Rcdmann 
(1976) found for a grassland canopy that the reflectance was considerably higher during the after- 
noons than in the morning for the same solar zenith angle. They propose that a portion of the 
observed variation is caused by canopy geometry in which preferential azimuthal orientation of , 
the leaves was assumed due to a preferential wind direction. The authors also suggested that plant 
water Content, leaf rolling, and atmospheric conditions may also be involved. If such variations in 
canopy geometry caused the above phenomena, the variations in canopy geometry due to temporal, 
spatial, stress, and wind variations may also account for the contrasting diurnal reflectance tra- 
jectories seen in this study and others. 

Variations in reflectance trajectories can be explained by other subtle geometric effects. For 
example, the simulated SRVC p£ results of an erectopliile canopy showed that for LAI values of 
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1.0 and 4.0 ^decreased with increasing solar zenith angle (Figure 3a). However, the opposite trend 
was observed lor an LAI of 7.0. The radiation transfer within vegetation canopies for highly absorp- 
tive wavelengths such as the simulated wavelength used in this study (0,68pm) is largely controlled 
by the probability of gap (PGAP) in the canopy as a function of source angle. Thus, the above 
variability can be explained as follows. At low LAI the PGAP to the ground varies greatly as a 
function of view angle. This variation for high LAI is small (Tabic 1). The spectral reflection (p^) 
and transmission (r^) for the canopy components and ground werep^= 0.8, t ^ = 0.4, and p^ = 

0.1 1, t ^ = 0.00, respectively (Kimcs, et al., 1979c). Thus, ground reflectance was significantly 
higher than the canopy components. The PGAP for a small zenith angle and a LAI of 1.0 was 
high. As a consequence, the observed p£ for small solar zenith angles should be high due to the 
relatively high contribution of ground reflectance. At large solar zenith angles the PGAP to the 
ground was small and the observed p£ is low due to the relatively small contribution of ground 
reflectance. For an LAI of 7.0 the PGAP was relatively constant with zenith angle and as a 
consequence the above effect was dominated by other geometric interactions. 

p 11 Reflectance 

Although p !l reflectances were not measured directly in this study, a large amount of data 
has been collected for vegetation canopies by other researchers. These reflectance measurements 
are compared with the results of this study. It is important to understand the inherent differences 
between p 11 and p c measurements. 

The simulated results within this study showed that the p^ and p^ 1 for the erectophile canopy 
as a function of increasing solar zenith angle decreased and increased, respectively. These differ- 
ences were explained by the anistrophic bi-directional reflectance function of the •simulated vege- 
tation canopy. Equations (1), (2), and (4) showed the p 11 and p c dependence on the bi-directional 
reflectance function. According to measurements from Coulson (1966) and FCriebel (1978) and 
simulated results from Oliver and Smith (1974), the nadir spectral reflectance was generally lower 
than all off nadir spectral reflectances for any given source direction. As a consequence, one 
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would expect the nadir p£ to be lower than the corresponding pj^ 1 reflectance. Such information 
is important for several applications. For example, accurate estimation of p 11 for global surfaces 
is important in calculating the planetary heat budget, iiaton and Oirmhirn (1979) have developed 
coefficients to estimate bi-hcmisphcrical measurements from hembphcrical-conical measurements 
for several target types. The purpose of these coefficients were to improve the prediction of p 11 
from the conical radiance values from space platforms. 

Coulson and Reynolds (1971) have measured p^J as a function of solar zenith angle of six 
vegetation canopies for six discrete wavelengths in the visible and near infrared regions. The 
authors found that the reflectances of most surfaces reached a maximum at solar zenith angles of 
80-70° and in general p” decreased with decreasing solar zenith angle, the authors suggested 
that these trends are a result of the effects of both canopy structure and the changing ratio of 
direct to diffuse light throughout the day. The magnitude of these independent effects on p" 
are not known. 

There are a large number of studies which showed a decreasing broad band bi-hemispherica! re- 
flectance (p 11 ) with decreasing solar zenith angle for vegetation canopies (Rijks, 1967; Davies and Butt- 
imor, 1969; Idso, ct al., 1969; Proctor, et al., 1972; Ripley and Rcdmann, 1976; among others). These 
studies used sensors sensitive to the visible and/or near IR regions. A typical parabolic curve from 
Proctor, ct al., (1972) is presented in Figure 5. As previously discussed, the changing p^ as a function 
of solar zenith angle was largely due to the non-lambcrtian nature of the bi-directional reflectance 
function, the changing solar irradiance field, and the directional field of view of the sensor. However, 
the parabolic relationship seen in Figure 5 for p 11 reflectance as a function of solar zenith angle cannot 
be explained by the above directional considerations since p 11 is a hemispherical measurement integrat- 
ing overall view angles. Other geometric effects must be operating. 

One explanation which is often presented is that the probability of a gap through the vege- 
tation components in the upper layers is generally lowest at the larger solar zenith angles. As a 
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consequence, less radiant flux will be attenuated and absorbed by the soil and lower layers of the 
canopy and a greater bMieinispherical reflectance will occur at the larger solar zenith angles. Tin’s 
is a simplified explanation of sonic very complex radiation transfers within vegetation canopies 
particularly in the near infrared region where strong multiple scattering occurs. 

Secondly, studies which measured broad band p 11 may introduce several source?) of error. 

For example, one source of error which is commonly not considered when measuring broad band 
reflectance values is the effects of variations in the spectral quality of solar irradiancc as a function 
of solar zenith angle and atmospheric conditions. When this variation is systematic in nature as a 
function of solar zenith angle, false reflectance trends can occur. As presented in the Appendix, 
simulated results showed that for a particular theoretical vegetation canopy the calculated bi- 
hemispherical reflectance (p !1 ) may increase as much as 25% from small to large solar zenith angles 
due to variations in the spectral quality of the solar irradiancc. 

Instrumentation error can also introduce systematic errors which cause false reflectance curves 
as a function of solar zenith angle. For example, using an Eppley pyranometer on an alfalfa 
canopy, Brown, et al., (1970) showed particularly that at greater solar zenith angles extraneous 
light (outside tire normal 2ir steradians) reached the sensor. The absolute reflectance error between 
a shaded sensor which reduced this extraneous light and an unshaded sensor is shown in Figure 6, 
This error can clearly introduce erroneous p 11 trends as a function of solar zenith angle. 

p u Reflectance 

In this article a specific p a function as presented by Kriebel (1976) lias been emphasized be- 
cause tile function was independent on the anistropic irradiancc field. However, there are a num- 
ber of studies which have made other types of bi-directional measurements on a variety of ma- 
terial types, for example, Egbert and Ulaby (1972), Eaton and Dirmliirm (1979), and Rao, et al., 
(1979) among others. 
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The trends and magnitudes of these studies suggested that these functions are highly variable 
and dependent on target types, geometric structure of targets, and the solar irrad lance field. These 
studies attempted to explain some of the variability in terms of the optical and geometric char- 
acteristics of the vegetation canopies. However, to date there has been a lack of general under- 
standing of the nature of the physical radiant interactions which take place in vegetation canopies 
in terms of the bi-directional function which is independent of the solar irradiance field, 

CONCLUSIONS AND RECOMMENDATIONS 

The wide variability seen in diurnal reflectance trends arc caused by variations in anistropic 
sky irradiance, canopy component geometry and optical properties, and type of reflectance meas- 
urement. In addition, systematic errors in in;* umentation and spectral integration can cause false 
diurnal reflectance trends. In most cases there is a lack of information concerning the effect of 
these variables on the magnitude of observed reflectance, 

Considering the large number of possible permutations of variables affecting diurnal reflectance 
trends, we believed that future studies should be designed to relate observed diurnal reflectance 
trends to the physical characteristics of the canopy and irradiance field. Collecting reflectance data 
only would be undesirable. Further decoupling the effects of the anistropic irradiance and the 
anistropic reflectance of the canopy on canopy reflectance measurements must be accomplished 
before a more complete understanding of diurnal reflectance trends can be obtained. 

It is unlikely that the variations In diurnal reflectances of vegetation canopies will be clearly 
understood until bi-directional measurements as suggested by Kriebel (1978) are more commonly 
derived and related to the geometric structure and optica! properties of the canopy constituents. 
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Table I 

Probability of gap (PGAP) Through a Theoretical I rectophile Canopy of 
Leaf Area Index (LAI) tqual to 1.0, 4.0, and 7.0 

Inclination View Angle Interval 

LAI 



0-10 

10-20 

20- 3C 

30-40 

40-50 

50-60 

60-70 

70-80 

80-90 

1.0 

0.0 

0.0 

0.0 

0.214 

0.392 

0.520 

0.613 

0.680 

0.720 

4.0 

0.0 

0.0 

0.0 

0.002 

0.024 

0.073 

0.142 

0.214 

0.268 

7.0 

0.0 

0.0 

0.0 

0.0 

0.001 

0.010 

0.033 

0.067 

0.100 
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l igure 3. SKVC simulated spectral hemispherical-conical rellectance p^ and bi- 
hemisphcncal rellectance p |* as a function of solar zenith angle for theoretical 
erectophile and planophile canopies. A photosynihetically active wavelength 
(O.t)K^m) was simulated for various leaf area indicies (I A I > . 
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Figure 4. Mean spectral hemispherical-conical reflectance as u function of 
solar zenith angle for three whrat stages lor Landsat bands 5 and 7. 
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Figure 7. Simulated bi-hemispherical reflectance from a theoretical 
vegetation canopy as a function of solar zenith angle. 
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APPENDIX 


APPENDIX 


When broad band bi-hemispherical reflectance (p H ) measurements arc taken as a function of 
view angle, error is introduced due to the varying spectral quality of the solar irradiance. When 
this variation is systematic in nature as a function of solar zenith angle, false diurnal reflectance 
trends can occur. The following simulation was performed as an example. 


First a basic distinction between broad band (e.g., entire solar spectrum) and discrete wave- 
length diurnal bi-hemispherical reflectance should be made. For very small discrete wavelength 
bands and ignoring path radiance and atmospheric transmission between the sensor and the target, 
the following ratio is measured: 

„ „ p” ‘ E \ ' R X < n 

• r * • a 


Where 

pH a bi-hemispherical spectral reflectance of target 
= hemispherical spectral irradiance 
= spectral responsivity of detector 
£2 = sensor related parameters 

However, when using a broad band (e.g., total solar spectrum) such as used by Proctor, et a!., 
1972 and Idso, et a.l., 1978, one is measuring; 
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Where; 

Xj ,X 2 = spectral limits of the effective detection of the sensor 

p 11 “ measured broad band bi-henrispherical reflectance for the (X t ,X 2 ) interval. 

For any given measurement period one wishes to measure a reflectance value for the entire 
spectral limits of the detector. However, it is clear that this measured value is dependent on the 
sensor’s spectral responsivity R\ and the solar irradiance E\. Indeed the F,^ is not a constant 
function but varies as a function of path length through the atmosphere and atmospheric con- 
ditions. Thus, it was hypothesized that diurnal reflectance trends as measured by broad solar 
band sensors may in some instances be only artifacts of the measurement procedure as reflected 
in Equation (5). 


The above hypothesis was examined in the following manner. Initially it was assumed that 
an ideal detector was being used with R\ = 1.0 units Xe(X!,X 2 ). In addition, for convenience the 
spectral solar irradiance was introduced as the product of the normalized spectral irradiance and 
the total irradiance. Thus, after the above transformations, Equation (S) becomes: 



Where: 

E n = normalized solar irradiance function 

* 


K 


dX = 1.0 


= total solar irradiance in the (X] , X 2 ) interval. 
The above equation can be further reduced. 


/ x 2 

CPS 


E?)dX 


( 6 ) 
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Theoretical E n functions were derived for several atmospheric path lengths as a function of solar 
zenith angle from data for a clear and dry atmosphere presented by Kondrat’yev (1965). The re- 
sulting curves are presented by Kimcs, et al,, (1979b). In addition, a typical reflectance curve of 
a leaf was utilized for p^ 1 as reported by Kimes, et al.; (1979b), This p” curve was held constant 
with changing solar zenith angle, thus isotropic reflectance of the canopy was assumed. Equation 6 
was then numerically integrated over 34 discrete wavelengths. The results for several solar zenith 
angles are presented in Figure 7. 

Note in Figure 7 that with a constant p” function we obtained an increase in total p" due 
to a changing E n function throughout the day. The shifts in the E n function throughout a theor- 
etical clear day weighted different portions of the p^ 1 function thus causing a changing total p 11 . 
Different responsivity functions (R^) of the detector could in some cases increase this effect. 

Such an artifact as discussed above can mask true reflectance changes of vegetational canopies. 
Gates, et al., (1965) have reported similar effects for individual leaf reflectance under clear and 
cloudy sky conditions. Although one may criticize the validity of the particular sky irradiance 
and plant reflectance curves utilized in this study, the above simulation does demonstrate that such 
effects can feasibly occur. The same effect can theoretically occur for broad band hemisphcrical- 
conical reflectance (p c ) measurements. 
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